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Abstract. We present spectroscopic observations obtained with the ESO Very Large Telecope (VLT) of seven candidate Lya 
emitting galaxies in the field of the radio quiet Q1205-30 at z=3.04 previously detected with deep narrow band imaging. Based 
on equivalent widths and limits on line ratios we confirm that all seven objects are Lya emitting galaxies. Deep images also 
obtained with the VLT in the B and I bands show that five of the seven galaxies have very faint continuum fluxes (I(AB)~26.8 
and B(AB)~27.3). The star formation rates of these seven galaxies estimated from the rest-frame UV continuum around 2000A, 
£f>^ , as probed by the I-band detections, as well as from the Lya luminosities, are 1-4 Mq yr -1 assuming a Hubble constant of 65 

I/-} ■ km s _1 Mpc -1 , fi m = 0.3, and Qa = 0.7. This is 1-3 orders of magnitude lower than for other known populations of star- 

forming galaxies at similar redshifts (the Lyman-Break galaxies and the sub-mm selected sources). The inferred density of the 
objects is high, 16±4 per arcmin 2 per unit redshift. This is consistent with the integrated luminosity function for Lyman-Break 
galaxies down to R=27 if the fraction of Lya emitting galaxies is ~70% at the faint end of the luminosity function. However, if 
this fraction is 20% as reported for the bright end of the luminosity function then the space density in this field is significantly 
larger (by a factor of 3.5) than expected from the luminosity function for Lyman-Break galaxies in the HDF-North. This would 
be an indication that at least some radio quiet QSOs at high redshift reside in overdense environments or that the faint end slope 
of the high redshift luminosity function has been underestimated. We find evidence that the faint Lya galaxies are essentially 
dust-free. These observations show that Lya emission is an efficient method by which to probe the faint end of the luminosity 
^ ' function at high redshifts. 
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al. 2000). So far it has not been considered efficient enough to 
seriously compete with the LBG technique, but it does have the 
advantage that spectroscopic confirmation is not limited by the 
broad band flux. Most recently the current success in identifi- 
cation of Gamma-Ray Bursts at high redshifts, has provided a 
completely independent selection technique for the host galax- 
ies of Gamma-Ray Bursts (e.g. Odewahn et al. 1998; Bloom 
et al. 1999; Holland and Hjorth 1999; Vreeswijk et al. 2000; 
Smette et al. 2001; Jensen et al. 2001; Fynbo et al. 2001). 

From the point of view of observations, what is still miss- 
ing is an understanding of the connection between high redshift 
galaxies selected in different ways. Fontana et al. (2000) find 
that photometric redshifts including IR colours select nearly a 
factor of 2 more high redshift galaxy candidates in the Hubble 
Deep Fields (HDFs) than pure LBG photometric selection 
does, to the same flux limits. However, this discrepancy can- 
not be resolved spectroscopically. While LBGs are selected by 
continuum flux, DLAs are selected by gas cross-section. Under 
the assumption of a scaling relation between the gas disc size 
and the luminosity for high redshift galaxies, and by normalis- 
ing this relation using the few observed impact parameters for 
z«3 DLAs, DLAs are predicted to be much fainter than the 
SLBG limit (Fynbo et al. 1999; Haehnelt et al. 2000; see also 
Ellison et al. 2001). Members of the population of galaxies pro- 
ducing DLAs are obviously not only found close to QSO lines 
of sight so we expect an abundant population of galaxies below 
the current spectroscopic limit for LBGs. 

From a theoretical point of view, the properties of this faint 
end of the high redshift LF is important in order to constrain the 
importance of stellar feedback processes (e.g. Efstathiou 2000; 
Thacker and Couchman 2000; Poli et al. 2001) and the fraction 
of the background of hydrogen ionising photons produced by 
high mass stars at high redshift (Steidel et al. 2001 ; Haehnelt et 
al. 2001). 

To probe the faint end of the LF at high redshifts we need 
methods to search for emission from high redshift galaxies 
fainter than R<25.5. Here several methods are possible: i) 
Using the Lyman break technique and photometric redshifts, 
LFs for z=3 galaxies has been presented by Adelberger and 
Steidel (2000) and by Poli et al. (2001) down to R=27. The 
faint end (R>25.5) of these LFs are uncertain due to the lack 
of spectroscopic confirmation, ii) Lya narrow-band imaging. 
in) Imaging of DLAs at very faint continuum levels with the 
Hubble Space Telescope (M0ller and Warren 1998; Kulkarni et 
al. 2000, 2001; Ellison et al. 2001; Warren et al. 2001; M0ller 
et al. 2001 in prep), and iv) deep searches for the host galaxies 
of well-localized (to within a fraction of an arcsec using the po- 
sitions of optical transients) Gamma-Ray Bursts. In this paper 
we focus on Lya emission as a method by which to study the 
faint end of the LF at z=3. 

In Fynbo et al. (2000, paper I) we reported on six faint can- 
didate Lya Emitting Galaxies detected in a very deep narrow 
band image of the z=3.036 radio quiet QSO Q1205-30 obtained 
with the NTT. Here we present photometric and spectroscopic 
follow-up observations of these 6 candidates (called S7-S12, 
detected at better than 5a) plus 2 marginal candidates (SI 3 
and S14, detected at the « 4a level). This paper is concerned 
mostly with the details of the observations, data reduction and 



the results concerning the physical properties of the Lya galax- 
ies. In a separate Letter we discussed the spatial distribution 
of the confirmed Lya galaxies and how it related to current 
models of structure formation in the early universe (M0ller and 
Fynbo 2001). In that paper we concluded that the Lya emitting 
objects are proto-galatic sub-units in the process of assembly, 
and for that reason chose to refer to them as "Lya Emitting 
Galaxy-building Objects" (LEGOs). Here, for consistency, we 
shall adopt the same acronym. The rest of the paper is orga- 
nized as follows: In Sect. 2. we describe the observations and 
the data reduction, in Sect. 3 we present our results, in Sect. 

4 we discuss our results in terms of star-formation rates and 
luminosity function compared to that of the LBGs and in Sect. 

5 we summarise our conclusions. Throughout this paper we 
adopt a Hubble constant of 65 km s _1 Mpc -1 and assume 
n m = 0.3 and fl A = 0.7. 

2. Observations and data reduction 

2. 1. Preimaging and construction of masks 

For the spectroscopy we used the Multi Object Spectroscopy 
(MOS) mode of FORS1 on the VLT-UT1 telescope. Pre- 
imaging for MOS was obtained in Service Mode during two 
nights in January 2000, well ahead of our Visitor Mode run, 
and consisted of 13 exposures each of 400sec in the Bessel B- 
filter. The 5 a detection limit of the final combined pre-image 
was B(AB)=26.7. Surprisingly this did not, however, allow the 
detection of continuum emission from all of the eight target 
objects S7-S14 (see Sect. 3.2 below). Furthermore, for one of 
those objects (S9) where we did detect continuum emission, 
there was a small but significant offset between the continuum 
position and the narrow-band position on the sky. For each of 
the eight target objects we therefore derived local transforma- 
tions between the NTT narrow-band image and the VLT B- 
band image, and placed artificial stars at the positions in the 
VLT image corresponding to the positions of S7-S14 in the 
NTT narrow-band image. Two of us did this independently and 
found positions that agreed to within 1/3 of a FORS 1 pixel, cor- 
responding to 0.07 arcsec. To minimize slitlosses, which could 
become severe if our invisible objects were not perfectly cen- 
tered in the slitlets, and/or if the objects were extended in Lya, 
we chose to use a slit-width of 1.2 arcsec for all slitlets. 

We chose grism G600B to cover the region of Lya at 
z=3.036 and grism G600R to look for other emission lines 
at longer wavelengths. Even though we did not expect to see 
any lines in the red spectra, it is imperative that they are deep 
enough that we can rule out the alternative interpretation that 
our objects are O II emitters rather than Lya emitters. G600B 
has a wavelength coverage from 3600A to 6000A (depending 
somewhat on the position on the CCD) and a spectral resolu- 
tion of 815 whereas G600R covers the range 5200A-7400A at 
a spectral resolution of 1230. 

We constructed three independent masks. In addidition to 
covering all eight candidate LEGOs, several of them in more 
than one mask, this also allowed us to obtain spectra of objects 
close to the quasar line of sight. For the mask construction we 
used the FORS Instrumental Mask Simulator (FIMS). We here- 
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Table 1. The log of VLT observations. 



date 


setup 


seeing 


Exposure time 






arcsec 


(sec) 


Imaging : 








2000 Jan 12, 17 


Bessel B 


0.68-1.02 


5200 


2000 Mar 5 


Bessel I 


0.55-0.75 


3750 


Spectroscopy : 








2000 Mai' 5 


MOS, maskBl 


0.78-0.96 


7200 


2000 Mar 4 


MOS, maskB2 


0.78-1.18 


7200 


2000 Mar 4 


MOS, maskB3 


0.61-0.66 


7200 


2000 Mar 4 


MOS, maskRl 


0.77-0.93 


5400 


2000 Mar 5 


MOS, maskR2 


0.57,0.78 


3600 


2000 Mar 5 


MOS, maskR3 


0.59-0.95 


5400 



after refer to the three G600B masks as maskBl, maskB2 and 
maskB3, and to the three G600R masks as maskRl, maskR2 
and maskR3. 

2.2. Spectroscopy and I -band imaging 

The spectroscopic observations were carried out during the two 
nights of March 4-5 2000 under photometric and good seeing 
conditions. We obtained total integration times of 4xl800sec 
for each of maskBl-maskB3 and 3x 1800sec, 2x3600sec and 
3 x 1800sec for maskRl, maskR2 and maskR3 respectively. As 
the seeing fwhm for all frames was significantly smaller than 
the slit-width the spectral resolution is determined by the see- 
ing and/or the size of the object. We also obtained 15x250sec 
images in the Bessel I-band filter. During the G600B observa- 
tions the CCD was binned 2 x 2 in order to reduce the influence 
of read-out noise. For all spectroscopic observations the CCD 
was read out in single port, high gain mode. For the imaging 
observations the CCD was read out using the four port, high 
gain mode. The journal of observations appears as Table [l]. 

2.3. Data reduction 

The pre-imaging B-band images and the I-band images were 
BIAS subtracted and flat-fielded using standard techniques and 
the individual reduced images were combined using the o- 
clipping algorithm described in M0ller and Warren (1993). The 
combined images reach 5(2)<r detection limits in 1 arcsec 2 cir- 
cular apertures of 25.9(26.9) in the I-band and 26.7(27.7) in 
the B-band (both in the AB-system). 

The spectroscopic frames were first BIAS subtracted. The 
flat-fielding was thereafter done in the following way. First we 
median filtered the flat-fields along the dispersion direction (x- 
axis) using a 30 x 1 boxcar filter for the G600B flat-frames and 
a 60 x 1 boxcar filter for the G600R flat-frames. Then we nor- 
malised the flat-frames by dividing the unfiltered flat-frames 
with the filtered flat-frames. Finally, we divided the science 
frames with these normalised flat-frames. The individual BIAS 
subtracted and flat-fielded science spectra were subsequently 
sky-subtracted in the following way. First we removed cos- 
mic ray hits from the sky-region of the 2-dimensional spec- 



trum by cr-clipping along each spatial direction column. We 
determined a 2-dimensional sky-frame using the background 
task in the kpnoslit package in IRAF. This sky-frame was then 
subtracted from the original BIAS subtracted and flat-fielded 
2-dimensional science spectrum. The individual reduced and 
sky-subtracted science spectra were then combined using cr- 
clipping for rejection of cosmic ray hits. 

1 -dimensional spectra were extracted using the apall task. 
For objects without or with very faint continuum we used bright 
objects from neighbour slitlets to determine the trace of the 
spectra. 

The 1 -dimensional spectra were wavelength calibrated us- 
ing the dispcor task. The rms of the deviations from the fits to 
3. order chebychev polynomia were 0.3-0.5A for the G600B 
spectra and 0.05-0.1 A for the G600R spectra. The high rms 
values for the wavelength calibration of the G600B spectra is 
due to the binning of CCD for the G600B spectra which imply 
fewer pixels per resolution element. 

3. Results 

3. 1. Spectroscopy 

3.1 .1 . Emission line candidates and line fluxes 

Our candidate list consisted of six "certain" candidates (Lya 
S/N>5) named S7-S12 and two "possible" candidates (Lya 
S/N«4) named SI 3 and S14. The candidates were selected on 
the basis of excess flux in the narrow-band as compared to their 
B and I band magnitudes (see paper I). A Lya emission line 
was detected spectroscopically for all objects except S14. The 
non-confirmation of S14 is consistent with the new and deeper 
B-band image also obtained during this run (Section 3.2), as 
the detection of B-band continuum for S 14 lowered the signif- 
icanse of the narrow-band excess below 3a. 

In Fig. [l] we show extractions of the blue grism spectra 
of all the candidates for which we confirm the presence of an 
emission line within the transmission curve of the narrow fil- 
ter. Only for S7 and S9 did we detect a faint continuum in the 
spectra. 

With variable seeing and 1.2 arcsec slits it was not possi- 
ble to obtain a spectrophotometric calibration. Nevertheless we 
obtained a rough calibration as follows. For each mask we had 
placed several slits on stars for which we had accurate B and 
I magnitudes. For the G600B masks we selected a star with 
a flat spectrum in B, and used that as a local standard for the 
flux calibration of emission lines. The calibration was done for 
each mask individually, and the results were combined after- 
wards for those objects that were observed more than once. The 
scatter for the objects observed more than once was about two- 
three times larger than the observational errors, confirming that 
this calibration is dominated by slit-losses as we expected. The 
resulting line fluxes are given in Table ^ (flux(slit)). 

The flux calibration described above is only valid if all the 
emission line objects are point sources. For extended objects 
the slit-losses will be larger than for the standard star, and an 
additional aperture correction must be applied. Comparing to 
the imaging line fluxes found with a circular aperture of diam- 
eter 3.5 arcsec (paper I), we find a mean aperture correction of 
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Table 2. Properties of the confirmed LEGOs S7-S13. Magnitudes are measured in a 3.5 arcsec diameter aperture. Two Lya 
fluxes are given, one as measured through the slit (flux(slit)) and one corrected to the 3.5 arcsec aperture (flux(aper) see text for 
details). The errors on the linefluxes are propagated errors from photon statistics relevant for determining the significanse of the 
detection, they do not include calibration errors which are estimated to be 30% because of the large aperture correction. 



Source 


wavelength 


redshift 


fwhm 


velocity width 


B mag 


I mag 


flux(slit)xl0 1Y 


fiuxCape^xlO 1 ' 


w obs 




(A) 




(A) 


km s 






erg s _1 cm -2 


erg s _1 cm -2 


(A) 


S7 


4911.6 


3.0402 


7.4 


<350 


25.61±0.11 a 


24.28±0.09 a 


1.04±0.15 


1.58±0.23 


> 61 


S8 


4911.1 


3.0398 


6.5 


<240 


~27.3 b 


~26.8 t 


1.70±0.18 


2.58±0.27 


456 b 


S9 


4905.2 


3.0350 


7.2 


<340 


25.36±0.07 


24.61±0.13 


3.51±0.18 


5.34±0.27 


> 164 


S10 


4905.6 


3.0353 


7.2 


<240 


~27.3 b 


~26.$ b 


2.09±0.26 


3.17±0.40 


456 6 


Sll 


4900.6 


3.0312 


7.2 


<440 


~27.3 6 


~26.8 6 


1.48±0.26 


2.25±0.40 


456 6 


S12 


4903.2 


3.0333 


8.2 


<520 


~27.3 b 


~26.8 6 


1.54±0.31 


2.34±0.48 


456 6 


S13 


4890.4 


3.0228 


7.3 


<450 


~27.3 b 


~26.8 6 


1.39±0.19 


2.11±0.29 


456 6 



a Including the nearby red object to the West of S7 proper 
b Measured in the median image of S8, S10-S13 




4750 4825 4900 4975 5050 4750 4825 4900 4975 5050 4750 4825 4900 4975 5050 

Wavelength (A) Wavelength (A) Wavelength (A) 



Fig. 1. The spectral regions around the Lya emission lines for the LEGOs S7-S13. For S8, S9, and S13 we show the spectra from 
several masks. As seen, the presence of emission lines is confirmed for all 7 candidates. For S14 (not shown) we did not detect 
an emission-line within the range of the narrow filter. 



0.46 mag. Line fluxes including this mean aperture correction 
are also given in Table || (flux(aper)). 

3.1 .2. Identifications and redshifts 

Before we can conclude that the emission lines are indeed due 
to Lya at z«3.04, we must first exclude any possible low red- 
shift "contamination". In the deep narrow band imaging dis- 
cussed by Kudritzki et al. (2000) 9 of 10 candidates were con- 
firmed as Lya emitters at z=3.1, while the last proved to be a 
galaxy at z=0.35 where the Oil 3727A line fell in the same 
narrow-band filter. There are other, albeit less likely, possible 
identifications (Mg II and Ne III). Those can all be tested by 
searching for additional emission lines (O II, H/3 and especially 



O III). For most of the possible interpretations the O III would 
be redshifted out of the G600B spectra. For this reason we also 
obtained spectra of all objects with the G600R grism. 

For none of the confirmed emission-line objects S7-S13 
did we detect any other emission line, neither in the G600B 
nor in the G600R spectra. In Fig. ^| we show the regions of the 
G600R spectra where the H(3 and O III lines would have fallen 
if the emission lines seen in Fig. |l] had been On at z=0.313 
(S7 at the bottom and S13 at the top). The spectra have been 
smoothed by a 7 pixel (7.4A) boxcar filter and regions with 
large errors due to strong sky-lines were set to zero. As seen, 
H/3 and O III lines are not present in any of the spectra. The 
question of possible alternatives to the Lya identification is an 
important one, and we shall return to a complete discussion in 
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Fig. 2. Shown here are parts of the G600R spectra were H/3 
and O III lines would have been if the emission line detected 
with the narrow-band imaging (paper I) were due to O II at 
z=0.31. S7 is at the bottom and S13 at the top. The spectra 
have been smoothed with a seven pixel (7.4A) box car filter. 
For S8, S9 and S13 the spectra are the average of 2, 3, and 2 
spectra respectively. 



Section 1.1 



, For now we shall assume that the lines are indeed 
due to Lya. 

The wavelengths and widths of the emission lines were de- 
termined by fitting them with Gaussian profiles. The uncer- 
tainty in the determination of the wavelength centroid is about 
0.3A. The results are given in Table || where we also list the 
resulting redshifts under the assumption that the lines are due 
to Lya. 

In order to constrain the intrinsic width of the lines we must 
know the spectroscopic resolution. An upper limit to the fwhm 
of the resolution profile along the dispersion direction can be 
obtained from the width of the slitlets and the dispersion. With 
a slit width of 1.2 arcsec (3 pixels) and dispersion of 2.4A per 
pixel the resolution would have been about 7. 2 A if the seeing 
had been worse than 1.2 arcsec. However, as the seeing was in 
all cases significantly better than 1 .2 arcsec the spectroscopic 
resolution is smaller than 7. 2 A fwhm. An upper limit to the 
spectroscopic resolution can be derived from the spatial pro- 
file of the spectra of point sources at wavelengths near 4900A. 
Converted to A the spatial widths are 5.9A, 5.7A and 4.8A for 
maskBl, maskB2 and maskB3 respectively. As the instrumen- 
tal resolution for FORS is slightly lower along the dispersion 
direction than along the spatial direction (T. Szeifert, private 
communication) these values can be used as lower limits on the 
spectroscopic resolution. From this lower limit we can obtain 
upper limits to the intrinsic widths of the lines by deriving the 



intrinsic widths that convolved with a Gaussian with a width 
corresponding to the lower limit on the spectroscopic resolu- 
tion reproduces the observed line widths. These upper limits 
are also given in Table ^. The line widths are smaller than 6A 
or 370 km s _1 for Lya at z=3.04. 

3.2. Imaging 

Fig. U and the two left panels in Fig. |] show regions of size 
12 x 12 arcsec 2 centred on each of the objects S7-S12 and SI 3- 
S14 from the combined NTT narrow-band image (top row), 
VLT B-band image (middle row) and VLT I-band image (bot- 
tom row). Only two of the objects, S7 and S9, are detected 
above our 2a detection limits of I(AB)=26.9 and B(AB)=27.7 
(detection limits for 1 arcsec 2 circular apertures). Those two 
galaxies were already detected in the deep NTT images pre- 
sented in paper I. For S9 both the B-band and I-band emis- 
sion is centred 0.5 1 ±0.05 arcsec south of the narrow -band po- 
sition. In the broad band images of S7 we see two nearby ob- 
jects. The center of the Lya emission is found 0.21 ±0.1 1 arc- 
sec east of the easternmost of the two, which we identify as the 
likely host of the Lya emission. We cannot be certain that the 
nearby western object is unrelated, hence the aperture magni- 
tudes for S7 given in Table || includes both of the objects. The 
not-confirmed (< 5er) candidate S 14 is clearly detected in both 
the B and I bands. The B magnitude is bright enough to con- 
firm that the apparent narrow-band excess in the old data set 
was not significant. 

For the five remaining objects we detect no broad band flux 
above 2a. In order to constrain the limit on the broad band 
emission further we registered the sub-images of the objects 
S8, S10-S13 to the centroids of their Lya emission and coad- 
ded them. In the right panel of Fig. ^| we show the median of 
the coadded images. We now detect a faint object in both the 
I-band and the B-band. To be able to compare the fluxes we ap- 
ply aperture corrections determined from a bright point source 
and arrive at I(AB) = 26.8±0.3 and B(AB) = 27.3±0.2 for the 
3.5 arcsec diameter circular aperture. 

If this flux was due mainly to one or two of the objects, 
they would have been visible on the individual images. Hence, 
we conclude that the flux must be fairly evenly distributed on 
most or all of the five objects, and that they each must have 
roughly the magnitudes measured in the combined images. It is 
interesting to note that unlike S7 and S9, the continuum emis- 
sion in the coadded frames is centered on the same position as 
the Lya emission. This suggests that the relatively bright con- 
tinuum object identified as S9 and the two-component object 
identified as S7 are, at least partly, due to a chance alignment 
of unrelated objects. Hence, the broad-band magnitudes given 
for S7 and S9 should probably be regarded as upper limits on 
their brightness. 

M0ller and Warren (1998) reported that on HST images of 
three Lya emitting objects related to the DLA at z=2.8 1 in front 
of PKS0528-250, there was evidence that the Lya emission 
was significantly more extended than the continuum sources. 
We measured the FWHM of the Lya emission of S9 and of the 
stacked Lya object, and found an intrinsic size (after correcting 
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Fig. 3. Contour plots of the combined NTT narrow-band (top), VLT B-band (middle) and VLT I-band (bottom) images for each 
of the six candidate emission line galaxies S7-S12 detected with S/N> 5 in the NTT narrow-band image (paper I). The size 
of the individual fields is 12x12 arcsec 2 and the fields are orientated with east to the left and north up. The contour levels are 
logarithmic. 



for the seeing) of 0.98 and 0.65 arcsec (FWHM). This supports 
the result reported by M0ller and Warren (1998). 

4. Discussion 

4.1. Line identification 

We now return in more detail to the identification of the de- 
tected emission lines. For a discussion of ways to discriminate 
between Lya and lower redshift objects for single emission 
lines around 8000-9000A, corresponding to very high redshifts 
(z>5) for Lya, see Stern et al. (2000). Here we discuss possible 
contaminants for z=3 Lya emitters. 

Mg II and Ne III can both be excluded easily since in this 
case we would have detected the stronger O II line. In the fol- 
lowing we will discuss how to reject the possibility that the 
observed line is due to O II. There are two independent ways 
to check the likelyhood of the identification as Lya rather than 
O II. One is by its equivalent width, the other by upper limits 
on line flux ratios. We shall here apply both tests, and as ref- 
erence point for the low redshift identification we have chosen 
the large samle of emission line galaxies from the survey of 
Terlevich et al. (1991). 

4.1.1. Equivalent widths 

The expected rest equivalent width of unextinguished star for- 
mation induced Lya can be as high as 200-300A (Chariot and 
Fall 1993; Valls-Gabaud 1993), which at a redshift of 3 would 
result in an observed equivalent width of order 1000A. Even 
with a fair amount of absorption by dust, the observed equiva- 
lent width would remain large. 



The continuum of low redshift O II galaxies is, in contrast 
to that of the Lya selected galaxies at high redshift, for the most 
part easily detected. In the sample of Terlevich et al. (1991) we 
find a median O II rest equivalent width of 59 A, and the 95% 
quantile is 163 A. In Table || we list the observed equivalent 
widths (W bs) of our lines. As discussed above, the broad-band 
aperture magnitudes of the two objects S7 and S9 probably in- 
clude flux unrelated to the emission line objects, and hence the 
calculated equivalent widths are lower limits. For the remain- 
ing objects it was necessary to coadd the images before we 
were able to measure their broad-band flux. For consistency 
we list the equivalent width measured on the stacked images. 

For an Oil galaxy at z=0.31 we would expect a median 
equivalent width of 77A and none of seven objects should be 
above the 95% quantile of 214A. This is incompatible with the 
observed distribution. 

4.1 .2. Line intensity ratios 

In order to compare the intensity ratios, we first had to deter- 
mine the upper limit for the non-detection of lines in the pre- 
dicted region for O III and H/3. For the flux calibration of the 
G600R spectra we used the exact same method as described 
in Section [D~T| forG600B. To maximixe the S/N for the non- 
detection of H/3 and O III emission lines, we stacked all the 
spectra. In Fig 5(a) we show the stacked spectrum (middle 
panel), and we have marked the expected positions of various 
emission lines. We find 2a upper limits to the log of flux ra- 
tios as follows: log(/om//oii) < -0.92, log(/ H /3/7oii) < 
-0.92, and log(/ No in//on) < -1-09. Above the stacked 
spectrum we plot the spectrum of the quasar in order to 
show the positions of the Lya and Civ emission lines. No 
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Fig. 5. a: Middle panel shows the stacked G600B (left) and G600R (right) spectrum of S7-S13. Lower panel shows, for com- 
parison, the spectrum of a z=0.224 emission line galaxy redshifted to z=0.313 so that the Oil line falls at the wavelength of 
the observed emission line of S7-S13. The stacked spectrum of S7-S13 has none of the lines expected if S7-S13 had been 
foreground emission line galaxies. Upper panel shows the spectrum of Q1205-30 to indicate the position of the C IV emission 
expected for AGNs. A C IV emission line is not detected, b: The box to the lower left of each of the figures (bl,2,3,4) marks the 
2<t limits for the flux ratios log(/oiii//on)> l°g(/H/3//oiiX an d log(/Noiii//on) under the assumption that S7-S13 are z=0.313 
galaxies with O II in the narrow filter. The observed flux ratios from the sample of Terlevich et al. (1991) are plotted ( x 's) under 
two assumptions used for the conversion of equivalent widths to fluxes ("mean" for bl and b2, and "worst case" for b3 and b4 
respectively, see text for details). Under both assumptions the 2a limits firmly exclude the z=0.313 hypothetis. 



Civ line is detected in the stacked spectrum. The flux ratio 
(log(/civ//Lya) < —0.92) excludes a significant AGN con- 
tribution to the Lya emission. Below the stacked spectrum we 
show for comparison the spectrum of an O II galaxy observed 
by chance during the same run in one of the slits we placed on 
random objects in order to fill the masks. This galaxy was ob- 
served at a redshift of z(0 ll)=0.224, but we have shifted it to 
z=0.313 for easy comparison. 

In Terlevich et al. (1991) we only found data for the equiv- 
alent widths of their emission line sample. In order to convert 
those to intensity ratios of two lines, we need first to deter- 
mine the ratio of the local continuum flux under those two lines. 
In all the published spectra of the Terlevich sample we there- 
fore measured the ratio of the continuum flux at 3727-3870A 
(Oil and Nelll) and at 5007A (OIII). We then determined the 
mean of this ratio and the maximum of the ratio (which is the 
worst case), and used both of them to convert the equivalent 
width ratios to flux ratios. The results are plotted in Fig. 5(b). 
In Fig. 5(bl,b2) we plot line flux ratios calculated from the 



mean continuum slope. The box and arrows in the lower left 
corner marks the 2a upper limits of our non-detections. In Fig. 
5(b3,b4) we repeat the same plots, but here we have used the 
"worst case" continuum for the conversion to flux ratios. Even 
in this case our 2a upper limits have no overlap with the O II 
galaxy distribution. 

In summary of this section we have applied two indepen- 
dent methods (equivalent widths and line flux ratios) to test 
how well our sample of emission line objects would fit if in- 
terpreted as low redshift Oil galaxies. Both tests reject the in- 
terpretation, and we conclude that all seven objects presented 
here are Lya emitters at z«3.04. 

4.2. Star-formation rates 

In paper I we calculated star-formation rates (SFR) based on 
the Lya fluxes. Here, for comparison, we calculate the SFRs 
from the continuum fluxes. As detailed in Section 3.2 we have 



reasons to believe that the continuum of S7 and S9 may be 
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Fig. 4. 7vvo Ze/f panels: Contour plots of the combined NTT 
narrow-band (top), VLT B-band (middle) and VLT I-band 
(bottom) images for the two candidate emission line galaxies 
S13 and S14 detected with S/N < 5 in the NTT narrow-band 
image. The size of the individual fields is 12x12 arcsec 2 (as 
in Fig. |^) Right panel: Same as Fig. || and the two left panels, 
but for the median image of S8, S10-S13 and with different 
contour levels. Broad band emission is detected in the median 
images at the exact position of the stacked narrow-band object. 



boosted by neighbour objects, therefore we perform the calcu- 
lation for the remaining five objects only. 

The restframe UV continuum in the range 1500A-2800A 
can be used as a SFR estimator if one assumes that the star- 
formation is continuous over a time scale of more than 10 8 
years. Kennicutt (1998) provides the relation 

SFRiMQyr- 1 ) = 1.4 x 1(T 28 x L v , 

where L v is the luminosity in the 1500A-2800A range mea- 
sured in erg s _1 Hz -1 . The observed I-band corresponds 
to the rest-frame UV continuum around 2000A, which falls 
well within this range. We can hence use the I-band flux 
in the median image of S8, S10-S13 as a measure of their 
SFRs. To derive L v we used the definition of the AB magni- 
tude to derive the observed flux (F v = io-oax(i(ab)+48.6) ) 
and finally the luminosity distance in our assumed cosmology 
(d/ um =8.58x 10 28 cm) to derive L v : 



F v x 47rd 2 „ m /(l + z) 



= 1.5 x 10 28 ergs~ l Hz~ l , 

where the factor (1 + corrects for the fact that L u is a 
specific luminosity (not a bolometric luminosity). Using the 
relation of Kennicutt (1998) we find a SFR of 2.1 M yr" 1 . 
This is 1-2 orders of magnitude smaller than the SFRs derived 
for LBGs (Pettini et al. 1998) and 2-3 orders of magnitudes 
smaller than for the sub-mm selected sources (Ivison et al. 
2000). 



From the Lya fluxes we found SFRs for the same objects 
in the range 0.3-0.5 hr 2 Mq yr _1 for a A = cosmology. 
Converting to the cosmology used in this paper this corre- 
sponds to 1.6-2.6 Mq yr _1 , which is identical to what we now 
find from the I-band flux. This result is inconsistent with the 
presence of large amounts of dust in those objects, which is not 
too surprising as one would expect that a targeted search for 
Lya emitters would preferentially find objects with very little 
or no dust. This does however imply that a large number of 
small star-forming objects at high redshift have essentially no 
dust in them. 



4.3. Space density 

The density of Lya emitters derived from the seven confirmed 
objects is 16±4 per arcmin 2 per unit redshift (seven objects 
within the 27.6 arcmin 2 field of view of the EMMI instrument 
of NTT and within the Az=0.016 range of the narrow filter). 
In Table [| we compare this to the results from other recent 
searches for LEGOs at z«3. It is seen from Table || that our 
survey found a larger space density than any other survey, but 
also that we reach the faintest detection limit of them all. The 
known LBG overdensity studied by Steidel et al. (2000) has a 
similar space density, but to a three times brighter flux limit. 
Down to their flux limit we would only have detected two of 
our seven sources. It is hence not clear from this if the volume 
around Q1205-30 is overdense in LEGOs, or if we simply see 
the effect of observing to a lower limiting flux. To be able to 
address this question, we need to compare our results to the 
extrapolation of the LBG LF. 

Adelberger and Steidel (2000) present a LF for LBG se- 
lected galaxies which is based on the HDF-North for the faint 
(25<R<27) and ground based LBG surveys for the bright 
(R<25.5) end. Integrating this LF down to R=27 leads to a 
predicted density of objects of 0.017 Mpc~ 3 h 3 . The comoving 
volume probed by our survey is 5.8 x 10 2 Mpc 3 h -3 . Based on 
the LF of LBGs we therefore expect ten galaxies in the volume 
and we find seven. However, only ^20% of the Lyman-Break 
galaxies show Lya in emission (Steidel et al. 2000). This may 
be an underestimate if Lya and continuum emission in general 
have different spatial distributions due to different slit-losses. 

With the Lyman-Break technique the only way to probe the 
R(AB)>25.5 part of the LF is to use the Hubble Deep Fields. 
This makes it impossible at present to study the faint end LF 
of any significant volume with this technique, and it is there- 
fore very uncertain. Indeed, significant differences in the num- 
bers of high redshift objects in the HDF North and South fields 
based on photometric redshifts have been reported (Fontana et 
al. 2000). The faint end of the LF is, however, important to 
study because a significant fraction of the star-formation, and 
therefore also the background ionising photons, may originate 
there. Based on the LF of Adelberger and Steidel (2000) there 
are roughly equal amounts of total luminosity from galaxies 
with R<25.5, 25.5<R<27 and 27<R<30. Hence, even if this 
LF is correct, only about one third of the total star formation at 
z=3 is traced by the LBGs studied in the ground based samples. 
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Table 3. The properties of other recent searches for LEGOs at z«3. 



Survey 


z,Az 


Area 

□ ' 


5(7 limit 

xl0~ 17 ergs -1 citT 2 


N 


dN 
dCldz 
#/□' 


field 


Confirmed 


Ref 


PKS0528-250 


2.81,0.019 


27 


3.7 


3 


5.9±3.4 


QSO, DLA 


all 


(1,2) 


HDFN 


3.43, 0.063 


29 


3.0 


5 


2.7±1.2 


blank 


all, 2 AGN 


(3,4) 


SSA22 


3.43, 0.063 


30 


1.5 


7 


3.7±1.4 


blank 


all 


(3,4) 


BR0019-152 


3.43, 0.063 


16 




7 


6.9±2.6 


blank 


3/7 


(3,4) 


Virgo 


3.15,0.043 


50 


2.0 


9 


4.2±1.4 


blank 


all 


(5) 


SSA 22a 


3.09, 0.066 


78 


3.0° 


72 


14±1.6 


LBG spike 


12/72 


(6) 


Q 1205-30 


3.04, 0.016 


28 


1.1 


7 


16±4 


QSO 


all 


(7) 



a The true 5<r level for this study is fainter, but only objects brighter than this limit were included in the catalog (Pettini, private communication) 

References : (1) M0ller and Warren (1993); Warren and M0ller (1996); (3) Hu et al. (1998); (4) Cowie and Hu (1998); (5) Kudritzki et al. 
(2000); (6) Steidel et al. (2000); (7) Fynbo et al. (2000) and this paper. 



In case the faint end LF needs to be revised upwards it will be 
even less. 

Assuming for now that the 20% Lya fraction can be extrap- 
olated to the faint end of the LF, we find that the volume we 
have surveyed has a comoving density of faint LBGs which is 
3.5 times that predicted by the HDF-North LF. This could indi- 
cate that at least some radio quiet QSOs reside in overdense en- 
vironments. Another likely interpretation is, however, that the 
faint end of the LF has a larger fraction of Lya emitters, in our 
case it must be of order 70% to fit the HDF-North LBG LF. A 
physical explanation could be that the objects in the faint end 
of the LBG LF are less dusty that the bright LBGs. 

4.4. Redshift distribution 

The Lya emitters were selected in a deep narrow-band search. 
It now remains to be tested if the objects span the entire width 
of the narrow-band filter, or if they cluster within a wavelength 
range smaller than the filter width. The width of the filter re- 
sponse was 23 A, and the full width spanned by the seven ob- 
jects is 21.2A. Monte Carlo simulations, where we randomly 
distributed seven objects weighted by the filter transmission 
curve and measured the resulting mean and std.dev., show that 
the seven redshifts of S7-S 1 3 are consistent (to within 1 a) with 
being drawn from a random distribution. Therefore, the struc- 
ture we have found is most likely larger than the redshift span 
we have covered with our filter. 

5. Conclusions 

We have reported on spectroscopic observations of eight candi- 
date Lya emitting objects at z«3.04. All six "certain" (> 5cr) 
candidates were confirmed, and of the two "possible" (< 5<r) 
candidates one was confirmed. To assess the most likely iden- 
tification of the lines we performed two independent detailed 
tests based on a large samle of low redshift O II galaxies. Both 
tests indicate that the only likely identification is Lya. This 
conclusion is strengthened by the fact that the original narrow- 
band imaging was centered on a quasar at the same redshift. 
The seven Lya objects are hence likely associated with the 
same structure as the quasar. 

The narrow-band images of the objects, as well as the large 
slit-losses of Lya emission we need to correct for when cali- 



brating the line fluxes, indicate that the Lya emission originates 
in extended objects. This is similar to the results reported by 
M0ller and Warren (1998) on Lya emission related to a DLA 
at z=2.81, where it was found that the Lya emission was sig- 
nificantly more extended than the continuum sources. As these 
authors point out, this could cause a severe underestimate of 
the Lya equivalent widths measured on spectra of high redshift 
galaxies obtained through a slit. 

An analysis of the redshift distribution of the seven con- 
firmed Lya objects shows that it is consistent with a random 
distribution in the redshift interval selected via the narrow- 
band filter. 

Despite deep detection limits only S7 and S9 were detected 
directly in the combined B and I band images. For S8 and S 10- 
S13 we registered the broad band images using the positions of 
the Lya sources and determined the median image of the five. 
This procedure allowed the detection of broad band emission 
at the level of B(AB)«27.3 and I(AB)«26.8 respectively. This 
means that S8 and S10-S13 belong to the faint end of the LF 
at z=3. The derived space density of LEGOs is consistent with 
the integrated LF of LBGs down to R=27. However, only about 
20% of R<25.5 LBGs are Lya emitters (Steidel et al. 2000). 
Hence, either the fraction of Lya emitters at the faint end of 
the LF is significantly higher than 20% or the space density of 
galaxies in the field of Q1205-30 is higher than predicted by 
the LBG LF. The latter would indicate that at least some radio 
quiet QSOs at high redshift reside in overdense environments. 

It has long been argued that Lya based searches for high 
redshift galaxies were doomed to failure, because even a small 
amount of dust would quench the Lya emission due to the res- 
onant scattering of the Lya photons. For this reason it has been 
virtually impossible to obtain telescope time for Lya survey 
work. We note here in passing that the recent very succesful 
Lya survey by Kudritzki et al. (2000) was aimed at low redshift 
planetary nebulae, and the work that supplied our candidate list 
was aimed at the imaging of a quasar absorber. It is virtually 
certain that neither of those programmes had been granted tele- 
scope time if the aim had been to search for a sample of Lya 
emitters at z=3. 

In this paper we find two strong arguments against the pres- 
ence of significant amounts of dust in the objects in the faint 
end of the high redshift galaxy LF. Firstly we find very large 
equivalent widths, in the upper range of theoretical predictions. 
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Secondly we find that when we calculate the star-formation 
rate from the continuum flux and from the Lya flux indepen- 
dently, we obtain identical results. Both of those observations 
are inconsistent with a dust-rich environment. 

The faint continuum magnitudes detected for S8 and S10- 
S 1 3 prove that Lya emission is a powerful method by which to 
probe the faint end of the galaxy LF at z~3. The next step is 
now to obtain statistically significant samples, containing sev- 
eral hundred Lya selected galaxies, in order to further constrain 
the properties of the faint end of the luminosity function. 
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